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ABSTRACT 

A recent lattice calculation of the QCD running coupling is presented. The coupling is extracted 
from the force between two static quarks in the framework of the valence quark approximation. 
A value of the A-parameter for zero quark flavours is determined: = 0.630(38)^ct = 

293(18)t^^ MeV. The first error is statistical, the second stems from the overall scale uncertainty 
in the string tension a. Combining this value with results from full QCD lattice simulations, 
we end up with the estimate A^ = 129(8)l6o MeV or ajjg{mz) = Q.imtfi- 



"'^Talk presented at the XVI Workshop on High Energy Physics, Sept. 14-17 1993, Protvino, Russia. 
^E-mail: bali@wptsO.physik.uni-wuppertal.de 



1 Introduction 



QCD, the theory of strong interactions, contains 8 free parameters, namely the quark 
masses, a (CP-violating) vacuum angle which is consistent with zero, and the running 
coupling parameter, that can only be fixed by experiment. For an empirical test of QCD, 
the following four steps have to be performed: 

1. measurement of 8 independent observables (e.g. particle masses), 

2. fixing the free parameters by use of QCD, 

3. predicting a value for another observable, 

4. comparison with experiment. 

It turns out that the second step is highly non trivial. So, most of this talk will deal with 
this task. I will use experimental data on the mass of the p-meson m^, and the string 
tension o", obtained from the quarkonia spectra, in order to predict the running of the 
couphng from QCD, which can then be compared to experiment. Since the precision of 
the lattice determination of as can compete with the experimental status quo, the result 
can be used as an input for phenomenological calculations. 

There are different ways to define a running coupling. Either one can state a value 
for Q;s(/i) at some fixed scale ix (like the mass of the boson m^-), the momentum 

at which the coupling has a given strength, or a A-parameter from the functional 
dependence a(g) = /(g/A), which can be obtained in perturbation theory, for high 
momenta transfer q. Here, I take a low energy scale /x, and calculate as{N x /x) on the 
lattice with N being a large number, in order to make contact with energies at which 
perturbation theory becomes valid. It is in this region where one can translate one 
scheme into another. 

The only known method which allows to treat QCD problems 

• nonperturbatively, 

• quantitatively, 

• from first principles, and 

• with control over possible systematic errors, 

is Lattice Gauge Theory. Unfortunately, it is a stochastical (non-exact) method, and thus 
statistical errors are introduced. Also, computational resources are a serious limitation. 
Since speed of supercomputers evolves according to an exponential scaling law in time, 
the latter problem seems to vanish. In the following, I will give a brief methodological 
introduction into lattice gauge theory, and comment on the advantages and caveats. This 
is followed by a section about how to determine the physical spacing of a given lattice. 
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a section about improved perturbation theory on the lattice, and a section about our 
actual calculation of as from the interquark force, and how to relate the result obtained 
to schemes which are more commonly used in perturbative QCD. At the end an outlook 
is given. 



2 Lattice Gauge Theory 

The first step in introducing the lattice formulation of a gauge theory is a rotation 
from Minkowski to Euclidean space (t it). Because of this rotation, information 
about real time dynamics is lost but static quantities like the mass spectrum remain 
unharmed. Then, space-time is approximated by a four dimensional hypercubic lattice 
with spacing a and extent La. tt/o serves as the required ultraviolett cutoff. After 
this is done, calculations of quantities of interest on the lattice are performed. Since a 
stochastic (Monte-Carlo) method is used, one also describes this procedure with the term 
"measurements are taken" . Finally, the results are translated back to the continuum. For 
this step, various extrapolations become necessary. The limit a — > (that corresponds 
to the limit ^ in asymptotically free theories) can be calculated by use of ordinary 
perturbation theory in the coupling parameter g^, once the perturbative region is reached. 
The infinite volume limit La ^ oo can be reached by use of Finite Size Scaling methods 
which are well established in statistical mechanics, and the extrapolation to small quark 
masses is controlled by chiral perturbation theory. 

In the following, I will briefly describe some basic notations, used in lattice QCD. 
The Wick-rotated action Sqcd = Sg + Sf contains a purely gluonic part, Sg, which can 
be written as 

Sg^-\[ d'xF^^F;^ = E (l - ^ReIV{C/^.(n)}) -rO{a') (1) 



Sw 



with (3 = 2N/g'^, and the number of colours = 3. is the bare lattice coupling. 
The discretizcd version of the gluonic action, Sw, is the so-called Wilson action. The 
continuum four- vector x is mapped to the discrete vector n = x/a. Instead of non- 
compact gauge fields, Afj,{x), SU{3) matrices are introduced: 



U^{n) - V 



I rn{a+ii) 

cxp ig j dxA^{x) 

\ Jna 



(2) 



In terms of these, the gauge symmetry is easier to preserve since products of such 
"link-variables" along a closed path, e.g. the Wilson action, are automatically gauge 
invariant. The lattice version of the Maxwell field strength tensor, then, is the or- 
dered product of hnk variables around an elementary square (plaquette) U^^{n) = 
exp {iga^ (F^,^(an) -|- 0{a))). In the following it will also be referred to as Un- 

Expectation values of observables O are calculated in the path integral approach as 
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their mean value over all possible field configurations in the finite box 




In order to solve this high dimensional integral, independent field configurations Cj = 
{Uj^\n)} are generated according to the measure dP{Ci) = (jln,f_idUj^\n)^ Q-fsSw{c^)_ 
The expectation value is calculated on a "representative ensemble" 

1 ^ 

(0)= lim T7E0(C.) . (4) 

Af— »oo JV T 
1=1 

This gives rise to a statistical error, AO, which reduces as is increased: AO oc 1/\/N. 

In the present calculation, the valence quark (quenched) approximation {Sf = const.) 
is used, in which vacuum polarization effects, caused by sea quarks, are neglected. In 
the language of perturbation theory this means that fermionic loops are dropped while 
gluonic diagrams are taken to all orders. The reason for this approximation are limited 
computational resources. However, all mass ratios from real world experiments are re- 
produced within an error of 6% [T[]. So, if full QCD with its dynamical fermions describes 
nature, discrepancies between full and quenched results are likely to coincide within 10% 
or so. This claim is also supported by the success of simple models like the naive quark 
model. 



3 How to set the mass scale 

In quenched lattice QCD, the inverse bare coupling (3 is the only free external parameter. 
So, the value of P determines the masss scale (and lattice spacing a): One measures a 
physical mass M in lattice units, and fixes the lattice spacing by the corresponding 
experimental value niphys- a{P) = mphys/M{P). On repeating this procedure for various 
P one can determine the dependence of the spacing on the bare coupling a^g"^), or, 
by inverting this relation, extract a running coupling g'^{a). As soon as a behaviour 
according to perturbation theory is observed (asymptotic scaling), this running lattice 
coupling can be translated into any other renormalization scheme perturbatively. 

The present simulations have been performed at the P values j3 = 5.5, 5.6, 5.7, 5.8, 5.9, 
6.0, 6.2, 6.4 and 6.8. L = 16 was chosen for the first five values, and L = 32 for the latter 
four. Below (3 = 5.5 (a ~ 0.25 fm) agreement with strong coupling expectations is found. 
Between 5.5 and 5.7 (a ^ 0.18 fm), as the (largest) correlation length in lattice units 
increases up to a few lattice spacings, a transition from discrete hypercubic symmetry 

to continuum rotational symmetry (j2i \Ri\ — ^ observed on the scale of a 

few a. Between 5.7 and 6.0 (a ~ 0.1 fm) ratios of physical masses, extracted from the 
lattice, remain constant within about 10%. From 6.0 up to 6.8 (a ~ 0.03 fm), these 
ratios remain constant within statistical accuracy of contemporary lattice simulations. 
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This is referred to as scaling in contrast to asymptotic scaling that does not set in for 
(3 < (9(10) at least. 

In order to fix the physical mass scale, we use the string tension cr = K/a^ which 
is defined to be the slope of the potential between heavy quark sources V{R) at large 
gg-separations 

K = -lun^F(R) , = . (5) 

The obvious advantages in extracting the mass scale from the potential can be summa- 
rized as follows: 

• The measurement of V{R) consumes little computer time since no matrices have 
to be inverted for calculation of propagators. 

• Finite size effects are negligible for La > 1 fm ^. 

• No extrapolations in the quark masses are needed. 

• The string tension is the most accurately measured dimensionful quantity on the 
lattice @, §. 

In the following I will discuss, which physical scale corresponds to ^/a and how to define 
it reasonably on a lattice with finite extent. 

Apart from investigations of Regge trajectories which suggest values ^/o^ ~ 400-450 
MeV, potential related information can be gained from the Charmonium and Bottomium 
spectra by integration of the Schrodinger equation |], If. It is practically only the 
potential at distances r between about 0.2 fm and 1 fm that determines these quarkonia 
spectra, since the corresponding wavefunctions almost vanish for smaller or larger qq 
separations. This is a lucky incidence because for smaller r relativistic corrections to the 
Schrodinger equation would have to be taken into account. For larger separations (r > 
1 fm), differences between the quenched potential (without sea quarks) and the real QCD 
potential are expected due to creation of a qq pair from the vacuum (string breaking). 
The phenomenological potential is well described by the Cornell parametrization 

V{R) = Vo-^ + KR . (6) 

Before discovery of the T', typical fit parameters have been ^/a = 455 MeV, and e = 
0.25 0. After inclusion of the T mass splitting, these parameters moved rapidly to values 
like v^=412 MeV, e = 0.51 § or v^=427 MeV, e = 0.52 ||. A careful inspection shows 
that the e ~ 0.25 and e ~ 0.5 parametrizations of the potential differ only weakly within 
the accessible range (0.2 to 1 fm). This is related to the fact that the fit parameters are 
highly correlated. 

We have attempted to fit our lattice potential (Fig. |lD for a fixed physical cutoff 
Ra > 0.3 fmQ, and different (3 to the Cornell form. In Fig. ^ we display the scatter of 

•^For f3 — 5.5 and (3 = 5.6, Ra > 0.5 fm, Ra > 0.4 fm have been chosen, respectively. 
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Figure 1: The (quenched) QCD potential V{R), measured at (3 = 6.2, (3 = 6.4, and 
j3 = 6.8 as a function of r in physical units. The mass scale is extracted from the string 
tension. 

fit parameters among different bootstrap samples in the i^-e-plane where we observe a 
correlation, similar to the ambiguity in parametrizing the quarkonium potential. To cir- 
cumvent this problem, Sommer has invented a more stable scale from the potential [jlO 



that is also well-defined in the screened case of full QCD. The idea is to calculate the 
dimensionless quantity X{Rq) = —F{Rq)RI from the interquark force. Since the po- 
tential models are most sensitive to the spectrum at a distance of about 0.5 fm, a scale 
tq = Rqq ~ 0.5 fm is defined by imposing the condition X{Rq) = 1.65. Using the fit 
parameters of Ref. one ends up with = 384 MeV. Ref. p| yields r^^ = 386 MeV 
and Ref. |^ = 402 MeV. Another set of fit parameters, taken from Ref. leads to 
= 416 MeV. We conclude with the value = (400 ± 15) MeV. 

Deviating from Sommer's prescription, we fit the potential to the Cornell ansatz and 
calculate 

' - {^e 

from the fitted parameters. All fits yield effective Coulomb coefficients consistent with 
the value e = 0.296, determined from the j3 = 6.2 data that carries the smallest errors. 
For convenience, we rescale the value of Tq ^ by the factor ^/h65^^J^296 = 1.164 in order 
to obtain a string tension cr^ = LIGAtq^. A bootstrap scatter plot in the ii'-e-plane is 
displayed in Fig. It can be seen that this definition of a string tension is insensitive 
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Figure 2: The fit parameter e versus K^^'^ (Eq. and the redefined string tension, 
calculated from Rq^ (Eq. The distribution 0/ 1000 bootstrap samples is shown. The 
higher the density of the points, the larger is the corresponding probability. 

towards the value of the Coulomb coefficient e. The string tension results are summarized 
in Tab. |l[ A detailed description of the data evaluation and fitting procedures will be 
given elsewhere Q. The physical value that corresponds to the redefined string tension 
is 

ar = 465 ± 18 MeV . (8) 

While ratios of pure glue quantities like the 0++ or 2+"*" glueball masses in units 
of the string tension exhibit good scaling behaviour above P = 5.7 [|ll|], the situation 
concerning observables with valence quark content is less satisfactory at present values 
of (3. In order to take these effects into account we have extrapolated our string tension 
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0.2204( 17) 


6.2 


0.386369( 
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6.4 
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Table 1: The average plaquette action, and the redefined string tension in lattice units. 
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values to f3 = 5.93 and (3 = 6.17, using the iJ-scheme that will be described below, and 
calculated the ratio of the p mass, taken from Ref. over the string tension at these 
/3-values and at /3 = 5.7. After this, linear and quadratic fits on the a dependence of 
this ratio have been performed, yielding a continuum value of 2 > mpj ^ > 1.7. This 
corresponds to 385 MeV < y/a < 455 MeV. The difference between this range and the 
number from the potential analysis (Eq. ^) might be caused by both, scaling violations 
still present above (3 = 6.0, and quenching. If we take this uncertainty into account 
and allow for an additional 5% quenching error, we end up with the conservative scale 
estimate 

a, = A65ttoo MeV . (9) 

By full QCD simulations the huge systematic error can probably be reduced to about 
5% or even less. In the following, 0".^ will be abbreviated as a, and arO^ as K. 



4 Asymptotic scaling and redefined couplings 
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Figure 3: The ratio y = a^^'^/Ajjg in the two-loop approximation, calculated by use of 
the bare coupling g"^ = 6/(3. The error band denotes the value y = 1.59 ± 0.10, extracted 
from the running of the coupling Uqq from the interquark force. 

The QCD /5-function can be expanded perturbatively 

-i^ = -bo9'-b,g' + --- . (10) 
dm a 
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Higher order terms are regularization scheme and gauge dependent. The first two coef- 
ficients are given by: 

llN-2nf _34:N^ - {13N^ -3)nf 

° ~ 3(167r2) ' ^~ 3Ar(167r2)2 " " ^^^^ 

Hf denotes the number of quark fiavours, N the number of colours. Integrating this 
/3-function leads to 

aA. = /(/):=e-^/^^-^(W)"'^'''° , (12) 

where A/, is an integration constant, not predicted by the theory. For sufficiently small 
g'^, the ratio y/a/A^ = y/K/aAi = \/K/f{g'^) is expected to remain constant. This 
perturbative behaviour is called "asymptotic scaling" . As can be seen from Fig. ^ there 
are large deviations from the asymptotic expression for (3 < 7. On the other hand it 
is the jS region (3 > 6 where good scaling between physical quantities like mass ratios 
is observed. So it might be that just the bare coupling is an ill behaved expansion 



parameter [0. The situation might be improved if one chooses a more appropriate 



renormalisation scheme by use of a more "physical" coupling parameter. 

Several suggestions have been made in the past for "better" couplings. I will review 
three such choices which have been proven to be successful, at least in the (3 range 
accessible to present day computer simulations. All are based on the idea of mean 
field improvement and were stimulated by Parisi [|T3|. The Fermilab group proposed an 



improved coupling ^fpNAL 05 IHl which corresponds to the continuum modified minimal 



subtraction (MS) scheme [IG] in the limit a ^ 0. A perturbative calculation [17 shows 



11 N'^ -1 

-d + 0.008204A^ + 0.00278n/ + ■ ■ ■ , d = . (13) 



ghsif^) g'if^)' — ^ ' 8N 



The large correction coefficient d, caused by tadpole diagrams, leads to the relation 
^Jis = T^e~^/'^'">AL = 28.81Ai for N = 3, Uf = 0. On the lattice, the average plaquette 



action (Sn) = {I — Un) can be expanded in powers of g'^ [0, [T^, 

(Sa) = c,g^ + C2g' + c^g'' + ■ ■ ■ . (14) 

The fields A^(x) are expected to fiuctuate around zero. So the mean field action is 
expected to be zero. Since the coefficient ci, that causes the deviations from this expec- 
tation, turns out to equal the large number d, that has been responsible for the huge 
difference between the two A-parameters, the idea is to substitute d = Ci hj {Sn)/g'^. 
This leads to the redefined coupling 

^ - + 0.008204A^ + 0.00278r2j (15) 



^FNAL(7i-/a) g^i-rr/a) 

with ApNAL = Ajjg. In Ref. |]15[ a different scheme is introduced, in which a coupling gy 
is defined from logarithms of small Wilson loops. 



Another empirical concept is the ii^-scheme, proposed in Ref. |21]. The starting 



point is the observation that large deviations from asymptotic scaling exist as well as 
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many quantities like the average plaquette are not well described by perturbation theory. 
The idea, now, is that the contributions, polluting the perturbative expansion of the 
plaquette, might be similar to those spoiling asymptotic scaling. By defining a coupling 
from inverting Eq. one forces the plaquette to behave according to perturbation 
theory, and hopes to find asymptotic scaling for dimensionful quantities as well. The 
E coupling is defined by gl = {Sn)/ci with Ae = e^2/(2cif>o)A^ = 2.0756Al for = 3, 
Uf = 0. Another suggestion 0, ^ to define a by inverting Eq. |14| after the second 
order. Then Ae2 = Al. 
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Figure 4: The ratio y = a^^'^/Ajjg in the two-loop approximation, calculated by use 
of the improved couplings versus the lattice spacing. The error band denotes the value 
y = 1.59 ±0.10, extracted from the running of the coupling Uqq from the interquark force 
(circles). The curves are just drawn to guide the eye. 

In Fig. the effective A^, obtained by use of the two-loop /5-function is plotted 
in units of the string tension versus the lattice spacing a for the above three effective 
couplings. We find that the behaviour is improved in all cases, though the i?-schemes 
appear to be more stable than the FNAL-MS-scheme. The differences between the val- 
ues, obtained from the effective schemes reflect still present deviations from asymptotic 
scaling. One might attempt to extrapolate y/a / Ajj^ to a = 0. If the corrections to 
the two-loop formula Eq. ^ can be resummed into an effective third order perturbative 
coefficient one might attempt to flt the data to the parametrization 

. ^ .(0) (1 - ^ .(0) (i + ,„(^(o)^(„^)J . M«) - x& (16) 
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where y{0) and C are fit parameters. 

Another possibihty that works equally well is to argue that the deviations are of 
nonperturbative origin and might be parametrized by a polynomial in the lattice spacing 
a. For 5.7 < /9 < 6.8 the data is well described by a linear fit, although from Eq. |16| it is 
clear that the behaviour might change rapidly as one approaches the continuum limit. 
As can be imagined from Fig. ^, with the described variety of fitting opportunities the 
extrapolated ratio might be anywhere within the range 1.4 < ^/o^ / A-jjg{0) < 1.8. 



5 Direct "measurement" oi as 



So far we have inserted a coupling g'^ (or calculated an effective g"^ at a fixed P), and 
measured the corresponding scale aig"^)- The effective couplings are also calculated from 
the plaquette, i.e. on a characteristic length scale a. This might cause nonperturbative 
deviations from asymptotic scaling that depend on the chosen discretization of the action. 
In order to reduce these "discretization errors", one should determine a coupling from 
observables with a typical scale of a few lattice spacings, instead. So far there have been 
three such approaches. 



1. Extracting couplings from logarithms of small Wilson loops by tadpole-improved 
perturbation theory p3|. 



2. Introducing a volume dependent coupling from the field response to the bound- 
ary conditions and exploiting finite size matching techniques on small volumes . 



3. Calculating a from a "physical" quantity like the interquark force [26, H, 22 



All these methods have their specific merits and disadvantages. Nonetheless, the results 
are fairly consistent with each other. In the first approach, it is not a priory obvious. 



what effective momentum corresponds to a Wilson loop with a certain extent [|15l . Also, 
one still encounters deviations from asymptotic behaviour that have to be examined 
carefully by simulating at various /?. An obvious advantage of this method is that Wilson 
loops are amongst the most simple quantities that can be calculated on the lattice. The 
second approach is a bit awkward since the boundary conditions have to be fixed. Thus, 
the generated configurations can only be used for the single purpose of calculating the 
coupling. While in the other two methods, can be computed from one simulation. 
Since the coupling is defined in a technical way (that allows for an easier perturbative 
expansion around the fixed background), the coefficient, connecting this scheme to the 
continuum motivated MS scheme, is large, and two-loop perturbative behaviour is only 
found for rather high energies. At present, it is not clear how this method can be 
generalized to full QCD with its dynamical fermions. The advantage is that the energy 
scale can be doubled with much less effort than when using a single-lattice method like 
approach 3, which I am going to discuss in the following. 
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The intuitive picture behind the concept of a couphng is that of local interactions. 
These are only realized at large energies where perturbation theory becomes a valid 
approximation. In principle, any dimensionless quantity that depends on a mass scale 
can be called a coupling. In the case of the quenched approximation to QCD, only one 
such mass scale exists. The /3-function, i.e. the logarithmic derivative of the coupling 
with respect to the mass, is universal up to two loops perturbation theory. Just the 
renormalization point differs from scheme to scheme. The schemes can be translated 
into each other by perturbation theory at large energies (small couplings). In the lattice 
formulation it is possible to define various nonperturbative gauge invariant couplings 
which can only be connected with each other in this high energy region. So one is forced 
to use a small lattice spacing, in order to make the results applicable in the framework 
of perturbation theory. A volume (La)"^ > 1 fm^ is needed to minimize finite size effects 
on the gg-potential 0, |^. Since the computer memory of our machine (an 8K CM-2) 
is limited to 256 MBytes L < 32 has to be chosen. These two constraints (the physical 
and the technical) set the limit a^^ < 6 GeV. Only potential values from R > \/2 are 
included into the analysis of the interquark force, in order to minimize discretization 
errors. As a check whether these are eliminated, simulations at different spacings a are 
needed. Fortunately, clearly perturbative behaviour is found for distances as small as 1 
GeV~^, i.e. over a scale factor 3-4 within the accessible energy range. 



0.1 



: • (1,0,0) 

0.08 - * ■ (1,',0) 

: ^ (2,1,0) 

" 0.06 - ° (1.1.1) 

- ° (2,1,1) 

^ (2,2,1) 

0.04 - m \ ' ' J 



0.02 - 



A Jll □ 



-0.02 - 



-0.04 - 



1 2 3 '1 5 6 7 8 9 

R 



Figure 5: Difference between the lattice one gluon exchange and the continuum 1/R. 

In order to increase the energy resolution we have to use information, obtained from 
qq separations where effects of broken continuum 0(3) symmetry are still visible. The 
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corresponding values of the force can be corrected from theoretical knowledge about the 
nature of these violations. The one gluon exchange, which causes a continuum potential 
V^(R) oc —1/R, is replaced by the Fourier transform of the lattice gluon propagator 



1 

R 



1 

R 



4vrG,(R) = 4. r cos(Rp) _ 



In Fig. 1^, the difference between the two tree-level results for all qq directions, realized in 
our simulations, is displayed. As can be seen, from R = 4 onwards it becomes negligible, 
and it is largest for the on-axis (1,0,0) separations. The one-loop coefficient has been 
calculated for the on-axis values only |]T^ , and gives — apart from a renormalization of 
the bare coupling parameter 



V{R) oc 



1 

R 



\n{AviR)R) 



(18) 



For large lattices and R > 1, the deviations of Ay{R) from the continuum Ay are less 
than 5%, and for the off-axis directions, they are probably even smaller. 

Exploiting this knowledge, we will reconstruct the continuum force F{R) and check 
our procedure by comparing results from four different lattice spacings reaching from 
2 GeV to more than 6 GeV. Then, we will define the coupling from the interquark force 

9N 

a,,{R) = -j^^F{R)R^ . (19) 



Finally, we calculate An = 1.048Atjc |28] from the running of the coupling according to 



a,g{R) = ^ {bo In (RaAn)-^ + b,/bo Inln (i^aA^)"') ' , (20) 

and translate our result into the MS scheme: aj^{q) = dgqir) — 0.082a^^{r) + ■ ■ ■ with 
r = 1/q. It is this small two-loop coefficient (0.082), and the dependence on just a single 
parameter (A^j), which turn the interquark force F[R) into a more user-friendly quantity 
than the potential V{R) itself is. 



5.1 Reconstruction of the continuum force 



We start from a parametrization of our data on the gg-potential which has been invented 
by Michael ^ 



V{R) = Vo + KR-^ + ^- 6ViR) , 6ViR) = el 
R R 



1 

R 



1 

R 



(21) 



The self energy Vq, the string tension K, the Coulomb coefficient e, the lattice symmetry 
correction-coefficient /, and / are fit parameters. / turns out to be Z ~ 0.6 for (3 = 6.0, 
6.2, 6.4, and 6.8, while / increases from 0.04 aX (3 = 6.0 up to 0.10 at (3 = 6.8 as it is 
qualitatively expected from the running of the coupling. The values have been found 
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Figure 6: Comparison between deviations of the potential at (3 = 6.4 from the contin- 
uum symmetry 5V{TC)/el, and [1/R] — 1/R, which has been used to parametrize lattice 
artefacts. 

to be acceptable for the fit ranges R > V2 {(3 = 6.0, (3 = 6.2) and R > {(3 = 6.4, 
j3 = 6.8). Tlie quality of the fits (for the example (3 = 6.4) is visualized in Fig. ^ where 
the theoretical value [1/R] — 1/R (circles) is compared against 5V{Y{)/{el) (squares), 
as computed from the data. As a check of the stability of the parametrization we have 
additionally performed four parameter fits by constraining / to zero. The parameters 
e turned out to be smaller in these cases, while / ~ 0.6 was remarkably stable. Only 
the accepted fit ranges changed to i? > 2 (/? = 6.0, (3 = 6.2), and > ^5 (/? = 6.4, 
P = 6.8), respectively. 

The force can be calculated from the potential by 

p^-^^ ^ V{K,) - V{K,) ^ ^ n2d/Rf) (22) 

with 

R = R{l + 0{i2d/Rf)) , i? = ^l±^ , d=^^ . (23) 

The smaller (larger) d is, the larger the statistical (systematic) errors become. The 
systematic errors stem from approximating the derivative by a finite difference. As a 
compromise, we have calculated the force from all pairs of potential values with distances 
restricted to 2 > 2(i > 1. In addition, an assumption about the smoothness of the 



a: 
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-0.02 



-0.04 - 



13 



resulting force has been made. As an illustrative example, let us assume V{R) = Vq — 
e/R + KR (Eq. 1). Using the above prescription with Ri = R — 1 and i?2 = -R + 1, we 
obtain F{R) = —e/ [R"^ — 1) — K. Since the expression is correct up to (9(4/i?^) only, we 
can substitute F{R) by F(R) with R := - 1 = R(l- 1/{2R^) + ■■■)• This leads 
to the correct derivative F{R) = —e/R — K. Exactly the same idea is used to calculate 
R{R,l, f/e) from the potential parametrization Eq. |2l]. This was done for / = /gt, as 
well as for / = 0. For I a systematic error Al = 1/5 has been allowed. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 7: The reconstructed continuum force from four (3 -values is displayed, together 
with the curve F{R) = —^aqq{R)/R^ for the parameter A/j = 0.66(4)0""'^/^. Horizontal 
errorbars denote systematic errors from taking the discrete derivative, vertical errorbars 
are statistical. The perturbative region is enlarged. Points with open symbols have been 
omitted from the determination of Aj^. 

In Fig. 1^ we display the resulting continuum force for all our /3-values. The system- 
atic errors from computing the discrete derivatives are included as horizontal errorbars 
while the statistical errors are shown as vertical errorbars. For small R the systematic 
uncertainties dominate. The success of our method is illustrated by the nice scaling 
behaviour over the four /^-values down to unexpected small distances (1.4 lattice units!). 
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/3 


i?-range 


akR 




6.0 
6.2 
6.4 
6.8 


1.9-2.0 
1.9-2.7 
2.0-3.7 
2.0-7.0 


0.1440( 43) 
0.1043( 25) 
0.0780( 24) 
0.0458( 52) 


0.653( 20) 
0.659( 16) 
0.661( 21) 
0.656( 76) 



Table 2: The plateau, in which a constant Aji{R) has been observed, and the averaged 
Ar in lattice units and in units of the string tension, respectively. 



Let me finally emphasize, that the fit to Eq. ^ has just been performed to estimate 
the relative weight of the one-gluon-exchange / = 0.60 ± 0.12. / is varied to account 
for systematic errors, arising from deviations from a Coulombic behaviour, due to the 
logarithmic running of the coupling. No additional information from the fits has been 
used for the numerical differentiation of the potentials. 
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0.08 
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0.065 
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2.5 
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Figure 8: The plateau in Aji{R) for (3 = 6.4. 



5.2 Calculation of the running coupling 

All data points on the force with Ri = 1 are excluded from our analysis to avoid un- 
controlled lattice artefacts. For all /5 and R values Ar{R) is calculated by inverting the 
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relation 



HR) = (iniRaAniR)) - ^ In (-2 ln(i?aA«(i?))) j , (24) 



obtained by combining Eq. 19 with Eq. 20. This ratio is plotted against R. A plateau 



is identified for Ra < 1 GeV~^ in all cases. As an example we show the (3 = 6.4 data in 
Fig. ||. The smaller errorbars indicate the statistical errors. The systematic errors are 
included into the larger errorbars. Below (Ra)^^ ~ 1 GeV the Aji{R) values increase a 
bit, and below about 0.5 GeV the values decrease towards zero. Since the data on the 
force is correlated, and systematic and statistical errors cannot clearly be disentangled, 
we do not attempt to fit the force but average the A^(i?) values over the plateau. The 
results are shown in Tab. ^ At /5 = 6.4, where Aij(i?) remains constant over a change 
of a factor two in the energy scale, and at /? = 6.8, where this factor is about four, 
these plateaus can clearly be identified. The 6.0 and 6.2 data yield consistent but less 
compelling results. So, we do not attempt to use them in our final analysis. Averaging 
the f3 = 6.4 and f3 = 6.8 results, and allowing an additional 5% systematic error on A, 
mainly caused by the one- loop lattice artefacts (Eq. |^), leads to 

Ar = (0.660 ±0.40)v^ . (25) 

The corresponding fit curve is plotted in Fig. ^. As can be seen, the data is well described, 
even for data points that had been omitted from the A-determination, in order to exclude 
discretization errors. On the other end, the data is in qualitative agreement with the 
perturbative prediction down to 0.5 GeV, surprisingly close to the pole Aji. 

The above value can be converted into 

= (0.630 ± 0.38)v^ = (293 ± 18l^^)MeV , (26) 

using the scale estimate Eq. ^ The second error is the overall scale error from the 
experimental uncertainty in ^/a. The four Ajjg values, extracted from the force at the 
different /3-values, as well as our final, averaged result (dashed errorband) are displayed 
in Fig. ^. The values scale so well that one might wonder whether we have overestimated 
the systematic errors. This question cannot be answered unless one-loop perturbative 
results on the off-axis potential are available. In Fig. ^, a comparison of the integrated 
perturbative force {Aji = 0.66^/o^) with the original, uncorrected potential is made. As 
can be seen, the data is described quite well, even without any lattice artefact corrections. 
Apparently, the theoretical curve has not been biased through the analysis procedure. 



5.3 Inclusion of dynamical quark flavours 

How can we relate our result, ajjg for the case of zero fermion flavours, to the real world 
that includes sea quarks? It is shown that our method works in principle and can be 
applied to dynamical conflgurations. In the meantime, it is possible to estimate the 
impact of dynamical quark flavours on the coupling by exploiting other lattice results. 
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Figure 9: The uncorrected lattice potential in units of the string tension, together with 
the integrated perturbative force (solid curve) for Aji = 0.66a^^^. 

Pioneering studies of full QCD have revealed that the sea quarks lead to a decrease of 
the lattice spacing a at fixed f3. However, dimensionless ratios of physical observables 
are in qualitative agreement with quenched results. The shift in j3, necessary to match a 
quenched lattice spacing a to the unquenched a, is called /?-shift. In the following, this /3- 
shift will be utilized to predict the corresponding change in the running coupling. From 
this matching, a large uncertainty in the energy scale arises but, since the coupling is 
running logarithmically, this has only a limited impact on the coupling at high energies, 
e.g. at the mass mz- 

For our extrapolation of ajju to mz, we use results, obtained by the MTc collabora- 



tion 1 29 1, for four flavours of staggered quarks with a physical mass of about 35 MeV. In 







5^^' 


5^°' 


(4) 2 

S'fnal 


(0) ^ 
fl'FNAL 




5.15 
5.35 


5.74(9) 
5.90(3) 


0.4760(1) 
0.4470(1) 


0.443 
0.418 


2.0597(4) 
1.8910(4) 


1.79(3) 

1.676(9) 


8.5i^:^ 

12.5l°| 



Table 3: Matching between four flavour QCD and the quenched approximation by use of 
the interquark potential. The data in the first three columns is taken from Ref. In 
the last column, the matching scale has been calculated. 
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this reference, the /3-shift has been computed by matching the dynamical qq potential to 
the static potential. The results are collected in Tab. |^. In column two the correspond- 
ing quenched [3 values are displayed. In columns three and four the plaquette values are 
collected. In columns five and six has been calculated in the FNAL-MS'-scheme. In 
the last column, the matching momenta vr/a in units of the string tension from our anal- 
ysis are displayed. They are comfortably within the perturbative region. A matching 
between the four fermion running coupling and the quenched coupling is done, yielding 
= (0.44 ± 0.10)A^. The scale vr/a cancels from this ratio. The error consists 
of a 12% part from the matching of both couplings at different vr/a, and a further 10% 
from possible scaling violations^, which might behave differently in the case of dynamical 
flavours (compared to the quenched case). 

Real nature has non-degenerate quark flavours. Thus, we have to incorporate sys- 
tematic errors emanating from the use of niu = rrid = nig = = 35 MeV into our 
prediction. What is the effect of shifting up to 1 GeV? If we evolve the coupling from 
1 GeV down to 0.5 GeV, an energy, where perturbation theory at least qualitatively de- 
scribes the running of agg, by use of the four flavour /^-function, evolve it back with the 
three flavour /3-function, and match it to the four flavour formula again (at 1 GeV), we 
find a decrease of A^ by 10%. A shift into the opposite direction of approximately the 



same magnitude |30] is caused by the fact that the light quarks have been too heavy on 
the lattice. So, we do not intend to change our central value but increase the systematic 
error and end up with 

= 0.44;J°A(^ = (129 ± 8+^g)MeV . (27) 

This can be converted into a value for five flavours with an additional error from matching 
the couplings at 0.75mb < /i < 2mb with nib = 4.8 GeV. The final result for five active 
quark flavours is 

«A«(^^^60%) = 0.102 (28) 
or, converting the scale error into an error on aj^g, this value is 



«Ms(^^) = 0.10201?^^ . (29) 
= Ai^ 

MS MS 

a-Ms{mz) > 0.115 from the lattice. However, the effects of dynamical fermions, as well 



Even, if we assume the extreme case A-^ = A^-L, it is impossible to obtain values 

' MS M5' ^ 



as the exact influence of quark masses [31 , have to be investigated carefully in the future. 



6 Summary and Outlook 

With a medium scale supercomputer simulation of the valence quark approximation of 
QCD, a running of the strong coupling "constant", as predicted by perturbation theory, 
has been found in the region between 1 GeV and 5 GeV. Qualitative agreement with a 

''As can be seen from Fig. ^, y/a/Ajjg, computed in the FNAL-Af S'-scheme, still differs by about 
25% from its asymptotic value at a > 0.1 fm w 0.25cr^^/^. 
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Source of uncertainty 


present error 


future error 


solution 




3% 


< 3% 


smaller a, higher statistics 


A^^ = Ar{R) 


5% 


0% 


one-loop off-axis potential 


v^= ? 


15% 


<5% 


dynamical fermions 


MS MS 


25% 


<5% 


dynamical fermions 


2% 


2% 





Table 4: Error sources in the calculation of AjJ^, and how these uncertainties can be 
removed/reduced. 



two-loop running was observed down to energies of about 0.5 GeV. From this, the QCD 
A-parameter (or alternatively the coupling at a given high energy scale) can be related 
to low energy observables like the hadron spectrum. Since the renormalization scheme 
in which the coupling was "measured" , is very close to the MS'-scheme, it is likely that 
the latter scheme gives reliable results down to energies as low as 1 GeV. Alternatively, 
one can use the lattice to define a renormalization scheme in a nonperturbative gauge- 
invariant fashion, e.g. the coupling from the interquark force. 

Combining the result with phenomenological potential models or data from lat- 
tice simulations of full QCD, a real world value can be estimated. The systematic part 
of the error, mainly caused by neglecting quark loops in the simulation, is dominant. An 
upper limit for ajjg{mz) is found: ajjg{mz) < 0.115. In Tab. sources of the scale er- 
rors, and how they might be removed in the near future, are collected. It should be notet, 
however, that, despite all systematic uncertainties, the accuracy of lattice simulations 
can compete with present day experiments. 

The method can immediately be applied to full (lattice) QCD. The computational 
effort of the present investigation has been moderate, compared to most other lattice 
studies of quantities of direct experimental relevance. So, an analoguous computation 
of the qq potential with dynamical fermions will become feasible in the near future. 
It should be amongst the earliest experiments on the first generation of TeraFLOPS 
computers. It will be very interesting to study the dependence of the running of as on 
the opening of new quark flavours. The impact of the quark masses and the transition 
behaviour at flavour thresholds can also be investigated. 
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